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4.9 GeV, the proton-proton contribution is larger than neutron-proton, but 
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1 Introduction 



Electromagnetic probes are very well suited for studies of strongly interact- 
ing particles under conditions ranging from free space to nuclear matter at 
high energies and densities. They carry away virtually unaltered informa- 
tion about the reaction since they do not suffer strong-interaction rescatter- 
ing. Our understanding of their properties continues to improve with the 
aid of both theoretical and experimental efforts. For instance, recent high 
statistics experiments have been carried out at Lawrence Berkeley Labora- 
tory by the Dilepton Spectrometer (DLS) collaboration in an effort to more 
firmly establish the mechanism of electron-positron production in 1-5 GeV 
nucleon-nucleon collisions The result of bombarding liquid hydrogen 
and deuterium targets with 4.9 GeV protons yielded a pd/pp ratio for e+e~ 
production that was nearly mass independent and ~ 2. This suggests that 
either bremsstrahlung is not responsible for the pair production yields or, if 
it is, that pp is roughly the same as np. However, conventional wisdom said 
that the dipole of the np system should be much more significant than the 
quadrupole of the pp system. Many dilepton production calculations have 
been done on pn, pA and AA systems at these energies using as a starting as- 
sumption that np bremsstrahlung dominates over pp [0-0. Even a complete 
Feynman diagram calculation supports this idea 0. By improving nucleon- 
nucleon bremsstrahlung calculations we settle this question. We also gener- 
alize to some of the inelastic channels: those of one-, two-, and three-pion 
final states. The cross section to produce pions in nucleon-nucleon scattering 
at these energies is comparable to or even larger than elastic scattering. One 
would therefore like to know if a proper treatment of the many-body elec- 
trodynamics gives a significant contribution. The pions produced would not 
only carry away some fraction of the momentum (thereby increasing pro- 
ton accelerations and therefore radiation) but they would travel very fast 
since they are relatively massless and if charged, would themselves radiate 
significantly. These two complementary features could result in rather large 
contributions to the low-mass dilepton spectrum. 

Other possible mechanisms for pair production in these nucleon-nucleon 
collisions are hadronic decays, two-pion annihilation and perhaps pion-pion 
bremsstrahlung. Pion-nucleon bremsstrahlung would also contribute but pre- 
sumably at a lower level. Dalitz decay of vr'' s contribute in a region where 
experimental acceptance is low, so in this study we neglect them. Contri- 
butions from Dalitz decay of the rj and radiative decays of the low-lying 
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nucleoli resonances might be important depending upon their production 
cross sections. We assume two-step processes for the resonances, namely, 
NN NR NNY (where R might be a A or an iV*). Guided by 
experimental data for the cross sections and analytically continuing the ex- 
pressions for real photon production to expressions for pairs of soft leptons 
using relativistic kinematics, we estimate these contributions. Two-pion an- 
nihilation was studied by Kapusta and Lichard in an attempt to explain 
possible structure near 2m.„ in earlier DLS data 0. What they found was 
that when acceptance filtered, the results were not peaked near threshold but 
rather near the p mass. The contribution near the peak was consistent with 
measured dilepton spectra when pt integrated. Finally, since multiple pion 
production is likely at these energies, we estimate the contribution from tttt 
bremsstrahlung in 4.9 GeV pp scattering. 

The paper is organized in the following way. In Sec. ^ we present the 
soft-photon approximation with which we calculate and the particular contin- 
uation to finite invariant masses we use. This will include partial corrections 
to phase space which are absolutely necessary in soft-photon approximations. 
Then in Sec. ^we apply the formalism to nucleon-nucleon collisions with very 
accurate parametrizations for the angular descriptions of elastic scattering. 
Section ^ contains an application to one-, two-, and three-pion final state 
many-body bremsstrahlung. The complicated matrix elements are approx- 
imated by their resonance structure since these inelastic channels proceed 
primarily through formation and decay of A-isobars. Dalitz decay of the 
rj s and radiative decays of the low-lying nucleon resonances are estimated in 
Sec. 1^. Inclusion of a model calculation of tt+tt" annihilation and its contribu- 
tion to the e~^e~ mass spectrum is shown along with our tttt bremsstrahlung 
estimates. Also in Sec. ^, we show the sum of the bremsstrahlung, rj Dalitz 
and A radiative decays and two-pion annihilation, as compared with dilepton 
yields in proton-proton collisions from the DLS. Then in Sec. we predict 
data for true cross sections instead of yields and present theoretical cross sec- 
tions without any experimental (acceptance) limitations. Finally, in Sec. |^ 
we summarize our findings concluding that many-body bremsstrahlung is 
the largest contributor to low-mass dilepton production in nucleon-nucleon 
collisions at these energies. 
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2 Soft-Photon Approximation 



One very pleasant feature of calculating within a soft-photon approximation 
is that the complicated algebraic structure of the diagrams separates into 
two somewhat less complicated pieces. Even if the final state is comprised 
of n > 2 hadrons, the matrix element for soft-photon production simplifies 
to the purely hadronic counterpart times a multiplicative function describing 
the complicated electrodynamics of the reaction ||^. Then as long as care is 
taken to approximately account for neglecting the momentum of the photon 
in the energy- conserving delta function of phase space, the differential cross 
section for quasi-elastic scattering while at the same time producing a virtual 
photon (electron-positron pair) of invariant mass M and energy go is 



a 



dM 



d^p+d^p. Rn{s') 



where Rn is n-body phase space 0, s' = s + — 2y/sqo and 



(1) 



cris = 



1=1 



nr=i d'p^ 



(go k ■ ^1 



(2) 



The square-bracketed expression in Eq. is the elastic {n = 2) or inelastic 
(n > 2) hadronic cross section and ggje ■ is the dimensionless electromag- 
netic weighting representing a coherent sum of all radiation fields involved. 
For this general process a + b l + 2 + ...+n, depicted in Fig. |I|, the 
electromagnetic four-current is 



QaVa QbPb 



Pa-q Pb-q 



+ E 



Pj-q' 



(3) 



where the Q s are charges in units of the proton charge. In this soft-photon 
method of calculating it is reasonable to take an angular average over the 
unobserved internal photon's solid angle. The result of the electromagnetic 



weighting is [ffU 
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<il\e-J? = -iQl + Ql + tQ') 

1=1 

-2QMl-^a-m0aJb) 

n 

+ Y.'^QaQ^{l-Pa■m0aJ^) 
i=l 
n 

+ J2'^QbQ^{l-Pb■m0tJ^) 
i=l 

n 

- J2 2QiQ, (1 - A ■ 4)X(A, 4), (4) 

i<j 

where the function 

[u ■ V — u"^ — VT^][u ■ V — v"^ — y/TZ] 
[u ■ V — u"^ + y/lZ ][u ■ V — v"^ + Viz ] 
and the scalar 

n={l-u-v)^-{l-u^){l-v^). (6) 

The /3 s are particle velocities p/E. This function presented in Eqs. 
is ultimately responsible for determining the relative strengths of pp and 
np bremsstrahlung. For 2 — > 2 hadron reactions, this function can be well 
approximated by a much simpler and much used expression when the mo- 
mentum transfer is small, e.g. it is proportional to t in np scattering for 
t < 4:771% 0. But we do not make such an assumption because it does not 
hold generally, a point that we beheve is not fully appreciated. We show 
in Fig. this electromagnetic weighting for 4.9 GeV np and pp reactions as 
a function of center-of-mass scattering angle. As we see, pp is suppressed 
only at the poles and rises to some rather large value for intermediate angles. 
For scattering angles in the forward hemisphere pp is larger than np. This 
is not inconsistent with the assumption that has prevailed in discussions of 
bremsstrahlung about np dominating pp since the assumption was based on 
a non-relativistic argument. At 100 MeV kinetic energy, for instance, the 
maximum value of the pp weighting (at cos^ = 0) is ~ 4.0 x 10~^, whereas 
the np weighting at that angle is ~ 7.0 x 10~^. For brevity, we do not show 
results at other energies, but already at 1 GeV the np domination breaks 
down. In the next section we parametrize pp and np angular distributions 
guided by experimental data, weight them with the functions plotted in Fig. ^ 
by evaluating Eq. (P, and obtain the e+e~ invariant mass spectra. 
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3 Nucleon-Nucleon Bremsstrahlung 



The absolute strength of nucleon-nucleon bremsstrahlung depends, among 
other things, on the differential elastic cross section da/dt. By definition, this 
is symmetric for pp scattering at any energy, whereas, it is slightly asymmetric 
for np at energies of interest to this study. This asymmetry can suppress the 
np bremsstrahlung contribution to dilepton production by a factor 4 for 4.9 
GeV reactions as compared with calculations using a symmetric da/dt [|TT| , 
. Clearly what is needed in order to minimize uncertainty in normalization 
are very accurate parametrizations of the experimental data. Functional 
forms that do quite nicely for energies 1-5 GeV are presented here. First for 
pp at 4.9 GeV we use 

^ = exp{Q{d-do)e{{7r-eo)-e){A + B\9-7r/2P) 

+Q{eo - e){c - D\e\^) 

+ Qi9-in-9o))iC-D\n-9\')} (7) 
where 6 is the Heavyside-Step function, 

B 
C 

and 7 = 2.3, 6 = 1.6, Oq = n/Q, d(Tpp{0)/dQ = 80 mb/sr, dapp{eo)/dQ = 0.9 
mb/sr and dapp{7T /2)/dQ = 5.0 x 10"'^ mb/sr; and then for np scattering 
again at 4.9 GeV we use 

da. 



i0o)/dn \ 




dn 



^ - exp{e{eo-e){A + B\e-eo\^') 

+ ei9-9o)iC + D\9-eo\')}, (9) 



. , danp{0o) 
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C = log 





da. 



danp{9o) / dVt ) 



dn 



and 7 = 2.3, 5 = 2.5, do = 99 degrees, danp{0)/dQ = 6.27 x 10^ mb/sr, 
danpiOo) / dVt = 1.62 x 10"^ mb/sr and danpij^) / dVl = 16.2 mb/sr. 



Comparisons are made with np |jT^, |T^, |T5[ and [0, data in Figs 



and |]b. Upon integrating these expressions over the scattering angle, or 
equivalently over momentum transfers — (s — Am%) < t < 0, we obtain cross 
sections of 11.9 mb and 11.1 mb for pp and np, respectively. Recall, there is 
a factor 1/2 in the pp case to avoid double counting of the identical particles. 

In order that we might compare to published data on dilepton production, 
we include the DLS acceptance filter to calculate a yield. The results are then 
acceptance filtered but not acceptance corrected |l[], and they have units of 
cross section. When Eqs. (J^) and (|^) are used in Eq. (0), and in turn in 
Eq. (|l|) (with acceptance filter and weight included), the resulting spectra 
differ by a constant factor of ~ 1.5. That the spectra differ by a constant 
is a consequence of our soft photon approximation. We show the yields in 
Fig. H, where the comparison is absolute. Proton-proton bremsstrahlung is 
stronger than neutron-proton at this energy! If we construct the ratio of 
a{np + pp)/a{pp), we find it to be 1.67. Therefore, it is not the measured 
ratio of order 2 that eliminates bremsstrahlung from being responsible for 
the e~^e~ yields; one must compare the distributions separately. Upon doing 
so, one if forced to conclude that simple bremsstrahlung is not responsible 
for the observed pp yields. One naturally thinks to estimate the inelastic 
channels' contributions to e'^e~ production. After all, at 5 GeV the cross 
section is mostly inelastic. 



4 Multiple-Hadron Final States 

The cross sections to produce up to and including five pions have been mea- 



sured in 5.5 GeV/c proton-proton collisions |T8[. The one, two and perhaps 
three pion final states are important for our considerations here. Their cross 
sections are 10.8, 11.4, and 11.7 mb, respectively. At the same time the 
elastic channel's cross section is 11.3 mb. So clearly, these inelastic channels 
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cannot be ignored. What is needed are the NN — >■ NNir, NN — >■ NNirir, 
and A^A^ NNninr matrix elements. One could of course perform a field 
theory calculation to get a fairly complete description of them. The in- 
teractions of vr and p mesons coupling to a nucleon and a delta are well 
known. However, the A propagator presents some technical difficulty mak- 
ing a full-blown calculation extremely lengthy. Furthermore, even a detailed 
one-boson-exchange calculation is not guaranteed to produce correct angular 
behavior for the differential cross section. We pursue a different approach. 
The process pp — * np-K^ proceeds primarily through or excitation of 
the nucleon, with subsequent decay into a nucleon and the positively charged 
pion. The corresponding diagrams are shown in Figs. §a and ^3, with ex- 
change diagrams not explicitly drawn. To this order, there are also effects 
from A° excitation as in Fig. |^c. Even though the matrix element for these 
diagrams contains influences from force-mediating meson-exchange, the reso- 
nance structure of the A propagator alone accounts reasonably well for what 
is observed in terms of invariant mass distributions. Therefore, we write the 
matrix element for the diagrams in Fig. |^ as 



-'^pp— >rip7r+ 0^ 



1 {Pi + PsY - (rriN + 

y/2 {pi + PsV - "^A + ^"^ArA 
_ 7 {P2 + Psf - (rriN + m^Y 

3V2iP2+P3Y -ml + iruATA' 

where a + b — > 1 + 2 + 3 corresponds to particles p + p n + p + tt^ . 
With this matrix element, we are neglecting the neutral excitation diagram. 
We have verifled that its inclusion is not signiflcant to the flnal dilepton 
spectrum, at least within this approximation. Keeping complexity to an 
absolute minimum while capturing the essential physics, we choose to neglect 
contributions from diagrams like the one in Fig. |]c. Then the cross section 
can be calculated by 

2m% l-^l%4/ , d^pid^p2d?P3 
d<^= r- ^=77^HPa + Ph-pi-P2-p?,)-^ E^TE^- (12) 



This reduces to an integration over the four essential flnal state variables 
describing phase space. The normalization of the matrix element is adjusted 
in order to give a cross section of a{pp npn^) = 8.03 mb. Similar ap- 
proximations are done for the matrix elements of the processes pp —>■ ppn^, 
np — > nrni^ , np —>■ npn^ and np ppTi~ . The cross sections we use are: 
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2.77 mb for pp ppn^, 8.03 mb for np — nrnr^ and finally, 2.77 mb for 
both np npiT^ and np — ppn^ . Future model calculations are needed 
to improve on these matrix elements, but for first approximation they are 
sufficient. 

With these hadronic one-pion-production cross sections, we insert into 
Eq. (1^) and subsequently into Eq. (|l|) in order to calculate each channel's 
contribution to dilepton production. Three-body phase space can be written 
as 

i?3(s) = J ds2R2is,S2,m%)R2is2,m%,ml). (13) 

When Eq. (|I]) is evaluated for the n = 3 case under consideration and with 
the acceptance filter included upon integration over dilepton pt and rapid- 
ity, a yield is obtained that can be directly compared with the simple pp 
bremsstrahlung presented earlier. In Fig. ^ we show the two single-pion final 
state contributions superimposed with simple bremsstrahlung. The largest 
contributor is the channel pp — > npn^e'^e' . For low invariant masses, it is 
larger than simple bremsstrahlung by nearly an order of magnitude. This 
many-body contribution has been ignored in previous calculations for dilep- 
ton production because it was assumed small. It is clear that such an omis- 
sion is not justified. The same analysis applied to the np channels gives 
similar results. Charged nucleon reaction partners are of course important 
for bremsstrahlung. But we have discovered that since the inelastic cross 
sections are comparable in size to elastic, having a charged pion in the final 
state significantly boosts the radiation. 

We then proceed to estimate the contributions from the two-pion final 
states. There are four channels each for pp and np scattering. The same 
kind of approximations are made in order to arrive at the matrix elements. 
This time, however, there are two A excitations and the isospin factors are 
different. The general structure of the diagrams is shown in Fig. |^. Again, 
we are neglecting diagrams that contain delta excitations on the initial lines 
since they are not essential. The matrix element we use for the process in 
Fig. is 

(Pi + Ps)^ - {ruN + m^Y {p2 + Pif - {niM + m^f 
{Pi + PzY - + irnA^A {P2 + PiY - m% + im/sX/^, 
{P2 +P^f - {niN + m^Y 
{P2 + PzY - rriA + i'mA^A 



+ 
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{Pi +P4y - m\ + im^T/^, ' 



(14) 



where pi and p2 label the protons and ^3(^4) labels the positive (negative) 
pion. Similar approximations give the matrix elements for pp ppn^Tr^, 
pp npn^n^\ pp — > nmr'^Tr'^, and np — > npn^n, np npn^n'^, np — > 
pp7r~7r° and np — > nnn^i:^. For simplicity, we assume an isospin equivalence 
forcing a numerical value of 2.85 mb for the cross sections of all two-pion 
final states. Measured pp channels differ from this by only a small amount, 
so it is not unreasonable. 

The phase space correction this time is Ra{s') / Ri{s) , where 

-^4(5) ^ j j ds2dm2R2{s,S2,m2)R2{s2,rn%,ml)R2{m2,mlf,ml), (15) 
s- m.- 

with ^2 = (mjv + iTij,y, $2 = {^/s — mjv — ttIt^Y, m2 — {iriN + fnT^Y 
f^t = (aA ~ V^)^- Naturally, this correction has a tendency to suppress 
higher massed pairs. The suppression with two pions in the final state is 
greater than it is having just a single pion. Carrying out the necessary phase 
space integration we find the resulting contributions to e+e~ yields from two- 
pion final states to be intermediate between simple pp bremsstrahlung and 
the largest single-pion channel. There are enough channels, however, so that 
when added together the low-mass contribution is slightly larger than that 
from single-pion final state bremsstrahlung. 

Finally, there is the issue of m > 2 pion final states. Note that at 4.9 
GeV there is enough phase space to allow as many as 12 pions in the final 
state. The cross section to produce three pions in pp scattering at 5.5 GeV/c 
is 11.7 mb. By doing a proper treatment of the five-body electrodynamics 
and by approximating the matrix elements in our by now familiar fashion, 
we find the contribution important only to the lowest two or three invariant 
mass bins. Specifically, using a cross section of 2.34 mb for the channel 
pp npTT+TT+TT-, we find dn^^"' /dM ^ 8.0 x 10"^ fxb/GeV for M = 0.075 
GeV. Admittedly this seems relatively large; but the five-body phase space 
correction very strongly suppresses the distribution for increasing invariant 
mass, bringing it down to 4.0 x 10~^ //6/GeV at M = 0.500 GeV. There 
are five charge configurations for three-pion final states which must all be 
included. Assuming charge independence for the cross sections, we arrive at 
our aforementioned value of 2.34 mb for each of these channels. The relative 
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smallness of the cross sections for four- and five-pion production justifies 
ignoring them here. 

In Figure ^ we show the sum of zero-, one-, two-, and three-pion final- 
state hadronic bremsstrahlung. One should compare to data only in the 
low-mass region where this soft-photon approximation is best. Upon doing 
so, we may conclude that these inelastic channels are responsible for much of 
the low-mass dilepton yield. A previous dilepton calculation which included a 
comparison between the soft-photon approximation and a Feynman-diagram 
method in 4.9 GeV np scattering showed the soft-photon result smaller than 
the diagram calculation by a factor of 6 at invariant mass 300 MeV. The two 
results converged to the same curve (as they must) in the limit M — > 
So one is encouraged that bremsstrahlung might even account for the yields 
at low-to-intermediate masses. 

5 Other sources 

Hadronic decays are the next most likely candidate for a mechanism of sig- 
nificant dilepton production. Of the possible hadrons that might be excited 
or produced in nucleon-nucleon collisions at energies ~ 1 GeV, we find that 
the r] and the A contribute the most. The rho and omega are somewhat be- 
low the delta, so we ignore them here. The delta has a direct decay channel 
A N'-f with measured branching ratio 0.6%. The photon might not satisfy 
the Einstein condition = 0, but instead appear as a virtual (massive) pho- 
ton. The precise mathematical formulation of the analytical continuation 
we use is identical to the calculation published first by Gale and Kapusta 
in Ref. The A production cross section we use is the measured value 

of 5.0 mb. We have also checked the contributions from radiative decay of 
the first few low-lying A^* s and found them to be small compared to the A. 
Calculation of the eta proceeds in the same way as the estimate from Ref. [|I| 
with one exception. We use the upper limit for the eta production cross sec- 
tion: 0.5 mb. Then a differential cross section is assumed which is gaussian 
in laboratory rapidity and exponential in transverse mass. For more details 
see Ref. H. 

At 4.9 GeV, the cross section to produce pions is larger than elastic 
scattering. The 1, 2, 3, and 4 pion production cross sections in pp scattering 
are 10.8, 11.4, 11.7 and 1.6 mb, respectively. Since the reaction zone has finite 

^The acceptance filter we use is Version 1.6, a later release than the one used in Ref. |19|. 
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extent, there is clearly some nonzero probability for pion-pion rescattering. 
Furthermore, vector dominance would suggest that in such scattering the 
rho resonance accounts for most of the cross section. It is well known that 
the rho has an electromagnetic decay channel, so it is reasonable to estimate 
the contribution to dilepton production through two-pion annihilation. This 
problem has been solved by Kapusta and Lichard in a search for structure in 
the dilepton mass spectrum 0. They found a structure appearing at the rho 
mass (instead of near 2m.„) when the distribution was acceptance corrected. 
For details of this kinetic theory calculation see Ref. 0. 

Pion-pion rescattering with bremsstrahlung will also contribute to dilep- 
ton production. Although the contribution from vrvr — 7r7re"'"e~ might be 
rather large , we must fold in some probability for two (final-state) pions 
to scatter, and then sum over all possible scattering energies and charge con- 
figurations. This requires complete knowledge of the dynamics. However, an 
estimate can be obtained from the following relation 

fdaiy \ r . dP dai 



TVTT 



, dm , (16) 

dM . J dm \ dM , ' ^ ^ 

\ / rescattering 2m^+M brems. 

where the probability for a pion to rescatter off another (pion) with invariant 
mass squared = (p^n + ^712)^ in a 4.9 GeV pp collision is 



dm 




(17) 



The annihilation cross section o"^+^_ , including a vector dominance form fac- 
tor, is taken from Eq. (15) of Ref. [0. Finally, the kinetic theory calculation 
of Kapusta and Lichard gives an estimate of the differential annihilation cross 
section [[^, da'^+^^- / dm. Having the probability for rescattering, we take the 
pion-pion bremsstrahlung expression developed within a soft photon approx- 
imation in Ref. and evaluate the (invariant energy) integral in Eq. (p!6|) . 
In practice, this is a finite sum J2i ^^i{dP / dm) x ((icr/(iM)brems- The re- 
sulting contribution to dilepton production is comparable to simple proton- 
proton bremsstrahlung in the low mass region and intermediate masses while 
it drops rapidly for pair masses above ~ 600 MeV. Note that this represents 
a lower bound on the pion-pion bremsstrahlung contribution for two reasons. 
Firstly, we have neglected other charge configurations and secondly, the prob- 
ability calculated with Eq. {\L7\ j has experimental acceptance contaminating 
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the numerator. The contamination reduces the probabihty in the mass region 
near threshold and should, in principle, not be there. However, the effect is 
very small. Since the pion-pion result is comparable to simple proton-proton 
bremsstrahlung at this energy, it is clearly not dominant. On the other hand, 
it is not insignificant either and should be considered in future calculations. 
Without performing a detailed numerical simulation, this is as far as we will 
take such estimates. 

The dilepton spectra arising from (total) bremsstrahlung, hadronic de- 
cays and from two-pion annihilation are all presented in Fig. |^. Dalitz de- 
cay of the eta peaks at a mass ~ 325 MeV and at a level slightly above 
bremsstrahlung. Therefore, it is a crucial ingredient in the final spectrum. 
The delta result is bimodal since, like all our results, it is acceptance filtered. 
The two peaks appear at 200 MeV and at the rho mass. From this it seems 
the delta is less important than the eta since the peak at 200 MeV is well 
below bremsstrahlung, and the other peak is masked behind the contribution 
from two-pion annihilation. 



6 Cross Sections 

In the previous sections we have presented dilepton yields since cross sections 
have not yet been published. Now we include a short section in which we pre- 
dict the cross section data for 4.9 GeV proton-proton collisions, since they are 
forthcoming. The dominant processes, mechanisms and calculational meth- 
ods are exactly the same as we have discussed in the previous sections, with 
one exception. Here we perform an acceptance correction, and therefore the 



answer is a true cross section ||2T[. In Fig. [1^ we show the resulting invariant 
mass distributions of true cross sections for bremsstrahlung, hadronic decays 
and two-pion annihilation. The prediction is absolute. 

Another useful comparison is to present cross sections for the various 
mechanisms without any experimental limitations included, i.e. without any 
acceptance effects. This theoretical comparison no longer depends on a par- 
ticular experimental apparatus. In Fig. ^ we show such a comparison. The 
hadronically inelastic channels contribute significantly in the low-mass re- 
gion. 
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7 Summary 



We have shown by properly calculating the electromagnetic interference in 
np and pp scattering and by accurately parametrizing the differential cross 
sections, that the e"'"e~ invariant mass distribution from pp bremsstrahlung 
is larger than np at 4.9 GeV. The single, double, and triple pion produc- 
tion hadronic cross sections were calculated while approximating the ma- 
trix elements by their gross properties of A formation and decay. Knowing 
the cross sections for these inelastic channels in pp scattering, we calculated 
their contribution to the dilepton mass spectrum. The single-pion produc- 
tion channel (with a charged pion) contributed the most to low-mass pairs. 
Dalitz decay of the eta contributed roughly at the same level in the narrow 
window near its peak. Finally, the sum of contributions from simple and 
inelastic hadronic bremsstrahlung, Dalitz decay of the r], radiative decay of 
A, and from tt~^it~ annihilation, satisfactorily describes the measured pp dis- 
tribution. We conclude that bremsstrahlung is indeed the largest source of 
low-mass dileptons in these 4.9 GeV pp collisions, but that it comes from 
many-body bremsstrahlung. 

Although the agreement with experimental data is quite good, it is not 
perfect. We do not attempt in this paper to achieve a closer fit. Our main 
goal was to establish the importance of the channels with multi-particle final 
states in electromagnetic radiation calculations. We have neglected several 
factors that can be important for precise quantitative interpretation of the 
experimental data. We have no effects of form factors in our dilepton emis- 
sion many-body cross sections. In this sense, the curves shown in this work 
represent a lower bound only. However, we have also neglected interference 
effects in the radiative resonance decays. These should be carefully examined 
since they have been shown to be of some importance in both dilepton and 
photon 1^ calculations. Radiative decays from multiple A excitations have 
also been neglected here. Potentially important effects in the pion-nucleon 
channels have not been treated either. Calculation of these rescattering ef- 
fects requires complete knowledge of the kinematics which is outside the 
scope of this paper. 

The hadronic inelastic channels will now play a major role in the interpre- 
tation and understanding of experimental data at these energies. An immedi- 
ate consequence of our study is that high energy heavy-ion data will require 
the use of state-of-the-art many-body numerical simulations, where multi- 
particle final-state channels and their appropriate electromagnetic weighting 
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must be included. It is of great importance to study the rise of the contri- 
butions with the incident kinetic energy to carefully map out the threshold 
effects. We eagerly await the upcoming data from the DLS and from HADES, 
the European dilepton collaboration. While it is true that the proliferation 
of new channels will complicate the many-body problem, we can on the other 
hand state that it will no doubt contribute to its richness. 
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Figure captions 

1. Hadron-hadron collision with an n-particle final state where the arrows 
indicate momentum fiow. A coherent sum of radiation for all charged 
external lines is computed. 
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2. Dimensionless electromagnetic weighting as it depends on the center- 
of-mass scattering angle in 4.9 Gev nucleon-nucleon collisions. 



Parametrizations of np and pp differential cross sections in (a) and (b), 
respectively. Experimental data are from (a) Refs. |l^, and (b) 
Refs. O, PI irai. 



4. Dilepton yield from simple pp and np bremsstrahlung as compared with 
DLS data from Ref. for 4.9 GeV proton-proton (open circles) and 
proton-deuteron (solid squares) inclusive reactions. Our pp result is 
shown as the dashed histogram, np is shown as the dotted histogram, 
and an approximate pd = pp + np is shown as the solid histogram. 

5. Pion production diagrams proceed through formation and decay 
as in (a), through formation and decay as in (b) and through 
A° excitation as in (c). The exchanged meson might be a pion, a 
rho-meson, or some other boson. 

6. Dilepton yields obtained from A^A^ —>■ NNtt bremsstrahlung as com- 
pared with simple nucleon-nucleon bremsstrahlung, and with pp data. 



7. Two-pion-production diagrams proceed via A^A^ — > AA — > NNtttt. 

8. The sum of bremsstrahlung contributions from A^A^, NNn, NNtttt, 
and NNtttttt final states in pp scattering. Data are again pp results 
from the Bevalac. 

9. Dilepton yields from various mechanisms: bremsstrahlung is presented 
as solid squares, radiative decay of the A is presented as solid triangles, 
Dalitz decay of the t] is shown as open squares, two-pion annihilation 
is shown as open diamonds, and finally, the sum of all contributions 
shown as the solid histogram. 

10. Cross section (invariant mass distributions) from various mechanisms: 
bremsstrahlung is presented as solid squares, radiative decay of the A 
is presented as solid triangles, Dalitz decay of the t] is shown as open 
squares, two-pion annihilation is shown as open diamonds and finally, 
the sum of all contributions shown as the solid histogram. 
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11. Cross sections without any experimental acceptance included. Labels 

indicate the following mechanisms: (a) simple proton-proton bremsstrahlung 
in the solid line, (b)-(e) are ppir*^ in the dotted line, npir'^ in the short- 
dashed line, ppvT+TT" in the dot-dashed line and npn~^n+7i~ final-state 
bremsstrahlung in the long-dashed line; (f) is eta Dalitz decay and 
finally, radiative Delta decay is shown in curve (g). 
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